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As a generalization of the mass-flux based classical stream-tube, the concept of momen- 
tum and energy transport tubes is discussed as a flow visualization tool. These transport 
tubes have the property, respectively, that no fluxes of momentum or energy exist over 
their respective tube mantles. As an example application using data from large-eddy 
simulation, such tubes are visualized for the mean-flovif structure of turbulent fiow in 
large wind farms, in fully developed wind-turbine-array boundary layers. The three- 
dimensional organization of energy transport tubes changes considerably when turbine 
spacings are varied, enabling the visualization of the path taken by the kinetic energy 
flux that is ultimately available at any given turbine within the array. 



1. Introduction 

The notion of a stream-tube ( Batchelor|1967 Fay|1994 ) as a tool for flow analysis and 
visualization is particularly useful because it maintains constant volume or mass flux 
across sections. A popular example of stream-tube based dynamical analysis is the ac- 
tuator disk and ideal flow model of wind turbines, establishing the relationship between 
power extraction and fluxes of kinetic energy at the stream-tube inlet and outlet (see 
Burton et al. 2001). The classical analysis is valid for steady and ideal flow. However, 



in three-dimensional high- Reynolds number turbulent flows, the transport of mean mo- 
mentum and mean-flow kinetic energy is often dominated by Reynolds stresses. Then the 
averaged flux of mass being visualized by a stream tube is no longer representative of 
the transport of other flow properties such as momentum or energy. As a generalization 
of the classic stream-tube, we consider the concept of momentum and energy transport 
tubes, defined such that there is no average transport of momentum or mean-flow total 
mechanical energy over their respective tube mantles. These tubes enable us to visual- 
ize the trajectory of transported properties across the flow. We illustrate the approach 
for the case of turbulent flow through large wind farms, looking for an answer to the 
question: "Where does the kinetic energy come from that is ultimately extracted by a 
given wind turbine?" To this end, we employ data obtained from large-eddy simulations 



(LES) of fully developed wind-turbine-array boundary layers (following Calaf, Meneveau 
I&: Meyers |2010| , with different stream- wise and span- wise turbine spacings and wind 
turbine arrangements. 

In Section [2] the concept of transport tubes for mean-flow momentum and mean- 
flow mechanical energy is presented. Next, to illustrate the concept, some analytically 
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tractable laminar flow examples are provided in Section [3] Subsequently, in Section |4] 
we demonstrate the use of these tubes for the interpretation of momentum and energy 
fluxes in large wind farms. There, we first evaluate conventional averaged-flow stream 
tubes that pass through a target wind turbine disk and then investigate momentum and 
energy transport tubes for various turbine spacings in the wind farm. Further discussion 
is provided in Section [5] Conclusions are presented in Section |6] 



2. Mass, momentum, and energy tubes 

We focus on incompressible viscous, statistically stationary turbulent flows with con- 
stant density p. A classical streamline of a stationary mean flow field is commonly de- 
fined as a curve V parametrized by x{s) € F (s € M), for which x{s) x it = 0, with 
u = [ui,U2,M3] the mean velocity vector in a particular frame of reference. A stream 
tube is then constructed by selecting a closed curve C, which is not anywhere tan- 
gent to the velocity, and considering the bundle of all streamlines through that curve C 
( [Batchelor 1967[ ). Consider a volume of stream tube $7, bound by the tube mantle M, 
and two cross sections Ai and A2. The volume-integrated continuity equation leads to 
/ /^^ pUiTLi da; -I- J /^^ pUiTii dx — 0, where n is the outward directed normal to the stream- 
tube control volume. No mass flows through the tube mantle M since by construction 
UiUi = there. 

To construct momentum or energy transport tubes, which have the property that 
there is on average no exchange of momentum or energy through the corresponding 
tube's mantle, we consider the vector fields formed by the total flux of these quantities. 
The total flux includes advective, turbulent and viscous fluxes. For the transport of 
linear momentum, we consider a direction characterized by constant unit vector ^ and 
components Q i^s an example, C could be any one of the Cartesian unit vectors i, j 
or k). Hence, the ^-momentum is puiQ, and for statistically steady flow its transport 
equation is given by 



where /j represents the body force, and 



dp 
dxi 



(2.1) 



(2.2) 



is the flux vector fleld of linear momentum (per unit mass) in the ^ direction (the index 
"m" refers to momentum), v is the kinematic fluid viscosity {v = p/ p with p the dynamic 
viscosity), and Sij — {dui/ dxj+duj / dxi) /2 the mean rate-of-strain tensor. Also, a linear- 
momentum-transporting velocity fleld can be deflned according to: 



2vSij]C,i 



(2.3) 



A similar notion of a diffusion velocity has been introduced in the context of deterministic 
particle transport methods ( Hermeliiie||1989 ) , and later used in particle vortex methods 
for viscous flows ( jPegond fc Mustieles 1990 Grant & Marshall 2005 ). Some similarity also 
exists with the notion of Favre averaging for compressible flow, where a mass-transport 
velocity field is defined by dividing the mass flux by the average density |Favre et al.\ 
(|1976D;|Smits & Dussaugel (|2006D . 
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Constructing a tube based on F„i (or u„j), we now find 

pFmjrij dx + 1 1 pF^jUj da; = - / / / da; + / / / J^i da;. (2.4) 



A2 J J At J J Jn dxi jjJn 

No momentum is transported through the tube mantle, since by construction Fm iUi — 
there. As a result, on the tube's cross-sections Ai and A2, the flux of linear momentum 
is constant, except for integral effects of sources and sinks of momentum, Qdp/dxi, and 
/jCi in the tube. The pressure effects can also be written in terms of the pressure at the 
inlet, outlet and mantle using — JJJq d{pC,i)/dxi dx = J f^^ p Qn.i dx + J f^^ p Qni dx + 
IImP dic- We remark that a classical jet evolving at constant pressure may be 
considered a stream-wise momentum transport tube, since the momentum flux across 
its sections remains constant and no forces or momentum fluxes act on its mantle even 
though a mass flux crosses the mantle - see Section [3] for details, where an elaboration 
of momentum and energy tubes (the latter is defined below) is presented for some sim- 
ple canonical laminar-flow cases. Note that the momentum-flux vector and related tube 
geometry depend directly on the choice of the direction <^ in which linear momentum is 
defined. The tube geometry also depends upon the velocity of the reference frame. Like 
classic stream tubes, generalized transport tubes are not Galilean invariant. For illus- 
tration on wind-farm cases in Section [4j we focus on stream- wise momentum along the 
incident wind direction only. 

Similarly we consider mean-fiow energy-transport tubes, based on the transport equa- 
tion for mean-flow kinetic energy (pK = puiUi/2): 

Q Qub^p du^ 

^(p^^Kj) - " ^^"^"^ a^J ~ 2^iSijSij+uJ„ (2.5) 

where 



Fk.j = Kuj + [u'^u'^ - 2vS^J)u^, (2.6) 

is the total kinetic energy flux vector fleld per unit mass, and the kinetic energy transport 
velocity is u^ j — F^j/K. These vector flelds may be used to construct energy-transport 
tubes. 

A difficulty for the interpretation of momentum and energy tubes is the fact that 
the pressure gradient acts as a source term. This is less of an issue when the pressure 
gradient is only due to an external pressure difference, such that the gradients only 
relate to an external force, and power inserted in the system respectively. However, when 
local accelerations or decelerations impact the local pressure (e.g. near wind turbines, see 
below) , the interpretation of these sources is less natural. For the study of energy fluxes, 
this can be remedied by looking at mean-flow total mechanical energy tubes. To this 
end, we decompose the mean pressure gradient as Vp = Vpoo + Vp — — /oo + Vp, where 
may be an external driving force per unit volume. Denoting the total mean-flow 
mechanical energy per unit mass as 



E = UiU,/2+p/p, (2.7) 



its transport equation reads 



" ''"''"^ 1^ " '^pS.jS.^, + u,U, + /.^oo), (2.8) 
with the total mechanical energy transport vector field Fe defined according to 



Fe,j = Euj + «u' - 2vS^j)ui. (2.9) 
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Related to this, a transport- velocity field may also be defined as Fej/E. Constructing 
a tube based on F^j, we now find 

J J pFejUj dx + J J pFe^jUj dx = - JJJ (^pSijSij - pu[u'j^^^ dx 

U:{1, + h^oo) Ax. (2.10) 



n 

No mean total mechanical energy is transported through the tube mantle and the flux 
across sections of the tube is constant, except for sources/sinks of mean- flow kinetic 
energy by the distributed force (ui/J, by mean- flow viscous dissipation (sink), and due 
to production of turbulent kinetic energy —pu'^u'-dvu/dxj (typically also a sink of mean 
energy). For conservative force-fields, one of course also has the option of including it 
into the definition of E via its potential function. Examples of energy tubes for some 
simple laminar-flow cases are briefly discussed in next section. 

Finally, in the particular case of ideal (inviscid) and steady laminar flow, we have 
u[u'^ = and 2vSij = 0. Hence, it is obvious that u ~ = Fk/K = Fe/E, from 
which it follows that stream tubes, momentum transport tubes and energy transport 
tubes all collapse, as conventionally used in ideal-flow, stream-tube analysis. However, in 
turbulent flows Reynolds stresses can affect momentum and energy fluxes considerably, so 
that these different tubes may differ greatly. This is illustrated with applications to flow 
in wind-farm boundary layers in Section |4] Furthermore, we can remark that transport 
tubes for other quantities such as vorticity, helicity, temperature or elements of Reynolds 
stress may be derived accordingly. 



3. Transport tubes for some simple laminar flows 

To first illustrate the concept of momentum and energy tubes, they are briefly elabo- 
rated in the current section for a few simple canonical laminar-flow cases. 

3.1. Couette flow 

We first consider laminar Couette flow, with an along-boundary (horizontal) velocity 
profile given by u{y) — y{U/h) and vertical velocity v — 0. Then the momentum flux is 
given by the following two components: F^ i = uu^ = {yU/hY and -F'm,2 = uvjn{y) = 
—v{du/dy) ~ —v{U/K). Therefore the slope of the tangent lines of this vector field is 
given by dym/dx = Fm,2/Fm,i = Wm/wm = -{i'h/U)y:^ = -i?e~i(i/,„//i)"2 (where 
Re = Uh/v). Integration yields momentum lines of the form 



h 



3 xq 



-,1/3 



Re 



(3.1) 



Figure [T] shows the resulting shape of these lines for two values of x^/h — ±1 thus 
enclosing a (2D) momentum transport tube. We used Re = 10 and the lengths shown are 
in units of h. The dot-dashed lines show regular streamlines. The momentum transport 
lines can be interpreted as follows: the flux of x-direction (horizontal) momentum flux 
that enters aX x — ~h (A-B) is transferred through this tube downwards towards the 
solid wall [ym = 0), where it is equalled by the viscous drag acting between x = —1 and 
X = 1. By definition there is no net momentum flux of any type crossing the solid lines, 
and since the problem is steady and there are no further forces acting (e.g. pressure), 
the entire momentum flux is absorbed at the wall. Alternatively, one may regard (e.g.) 
the bottom horizontal (dot-dashed) line as the top wall that is being dragged from left 
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Figure 1: Couette flow: dot-dashed lines show sample streamlines, solid lines are momen- 
tum lines, while dashed lines denote kinetic energy lines. 



to right. Then the total drag force acting on this top wall between x = —5 and x = —3 
(segment A-C) is "transmitted" via the momentum transport tube towards the bottom 
wall as indicated by the solid lines. Near the wall, the transport velocity becomes vertical 
as more and more of the momentum transport occurs through viscous diffusion, while 
the momentum being transported vanishes. Because of the latter effect, the magnitude of 
the vertical transport velocity diverges to infinity, while the transport lines remain well- 
defined. We also remark that at increasing Re (or increasing h away from the bottom 
wall), the transport lines become more horizontal, as inertia in the horizontal direction 
dominates over viscous diffusion. Close to the bottom wall, viscous diffusion dominates. 



Next, we consider the kinetic energy transport flux, see Eq. (2.9). Similarly as with 



momentum transport, one can show that the energy lines are given by 

6 xo - 



f \ r « 1 1/3 



h 



Re 



(3.2) 



The resulting lines, starting out at the same points as the streamlines and momentum 
lines at a; = —5 in Figure [l] and shown using dashes, curve down towards the wall more 
quickly than the momentum lines. The entire flux of kinetic energy that enters the energy 
tube at X = — 5 (segment A-B) is transported towards the wall while being dissipated 
into heat inside the tube. Since no work is being done on the bottom wall, the entire 
energy is dissipated inside the tube before reaching the bottom wall. Conversely, there 
is work done by a moving top wall, e.g. along segment A-D, which is then is dissipated 
into heat inside the energy tube. 

3.2. Poiseuille flow 

Similarly, if we consider simple laminar Poisseuillc flow of the form u{y) = y{l — 
y)Gh? / {2v), with G — —{1/ p)(ip/dx and y non-dimensionalised with channel total 
height ft., we obtain Fm,i ~ u{y)'^ and Fm,2 = ^(1 — 2y)G'/i/(2i/). The momentum- 
hne slopes become dj/m/dx — — (8i?e?i)~^(l — 2y)[j/^(l — y)^]~^, where Reh = Uh/v and 
U = Gh?/{lQu) is the channel mean velocity. The slope depends upon the Reynolds 
number, with steeper slopes corresponding to lower Reynolds number (stronger diffu- 
sion transport) as expected. Integration yields momentum lines given hy y = ym{x) and 
passing through (a;o,yo)i according to: 

41n f^^V(22/~3)(2y+l)(l-2y)2-(22;o-3)(2yo+l)(l-2yo)' = ^{x-x^). (3.3) 
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Figure 2: Momentum lines for Poiseuille flow at three Reynolds numbers 




Figure 3: Representative streamlines (dot-dashed) and momentum lines (solid) for lami- 
nar round jet similarity solution. 

The results are shown in Fig. [2] for three Reynolds numbers. The interpretation is that 
the momentum added in the bulk of the flow through the pressure gradient is transported 
towards the side-walls as shown in the figure. We note that kinetic energy lines have the 
same shape but with twice the slope, and the work done by the pressure gradient is 
dissipated entirely, before reaching the walls. 

3.3. Laminar round jet and wake 
The similarity solution for laminar round jet is given by the stream-function ip — h'xf{ri), 
where rj — r/x (White 20061. The axial velocity is u = {v/r)f' and the radial velocity 



IS V = {v/r){rif' — /). The classical solution is fijf) = (cry)^ [l -I- {crj/2Y\ where c is 
related to the jet momentum flux and Reynolds number. The shape of the constant-i/' 
streamlines are visualized in Fig. |3] using dash-dotted lines. The horizontal momentum 



lines can be obtained using F„ 



I? and F„ 



uv 



vdu/dr = {vlrYf{^f' ~ f + 



1 — rjf / f"). The slope of momentum lines is then given by 

, f-ff-vf" 



drm 
dx 



F„ 



(3.4) 



since /' — //' — rjf" = for the round-jet similarity solution. As a result, the momentum 
lines are straight lines (solid lines in Fig. [s]) . This is of course expected, since momen- 
tum flux is constant in sections of cones bounded by a fixed 77, a basic requirement 
since f'^dr]' only depends upon rj. The lines in Fig. 3 are helpful in visualizing how 
momentum is being brought towards the outer entrained fluid. 

The kinetic energy lines in the far field (using the boundary-layer approximation that 
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K = ^u^ {v << u) also are lines but have twice the slope of the momentum lines (for the 
same value of rj). The "faster" spreading does not mean that comparatively more kinetic 
energy than momentum arrives at the entrained fluid, since part of the kinetic energy is 
dissipated. 

The far-wake solution with horizontal velocity u — Uq ~ Uc and centerline defect 
velocity of Uc — (C/x) exp(— ry^) with rj — r[Uo/{4:iyx)]^^^ leads to a radial inflow given 
hy V — {v /Uq)^^^C x~^/^ rjexp^—rj^). The radial momentum flux is i^m ,, — uv—udu/dr « 
Uqv + vduc/dr in the far-wake, which can be shown to lead to Vm — 2u. Similarly, the 
kinetic energy transport vertical velocity is given by vk = 3w. Hence, momentum and 
kinetic energy are being transported into the wake from the outside along steeper, but 
similar, influx "trajectories". In the wake, mass, linear momentum and kinetic energy 
are being replenished. 



4. Transport of mass, momentum, and energy, in large wind farms 

We now turn to the application of transport tubes for the visualization of momentum 
and energy transport in the flow through wind-farm boundary layers. 

With the increase in size of land-based and offshore wind farms the problem of farm 
performance is becoming an important research topic (Emeis & Frandsen'1993'; 'Frandsen 
~et a/..„2069l llvanell et aL ,2009; Barthelniie et aL||2010, Cal et a/.||2010, ,Lu fc Porte- Agcl 



2011 



Meyers fc Meneveau|2012 ). For very large systems, the notion of the asymptotically 



large (infinite) wind farm becomes relevant (Emeis fc Frandsen 1993 Frandsen et al. 2006 



Calaf et a/.|2010 1 . This limiting case can be conveniently studied in numerical simulations 



using periodic boundary conditions in the horizontal direction, as has been done in recent 
LES studies of wind farms ( |Calaf et a/.||2010[ [Meyers fc Meneveau||2012[ ). 

For a lone-standing turbine, physical mechanisms related to power extraction are rea- 
sonably well described using a stationary stream-tube analysis, neglecting effects of vis- 
cosity, and Reynolds stresses. Conservation of mass, the Bernoulli equation, and consid- 
ering differences in up-stream and down-stream momentum fluxes, lead to concepts such 
as the Betz limit for wind-turbine power extraction, wind-turbine momentum theory, etc. 
( Burton et aL|[2001 ). In real wind farms, however, wind-farm induced turbulence levels 
are much higher, so that turbulent fluxes become already as important as ideal terms 



in stream tubes that extend 2D upstream and downstream of turbines (Lebron et al. 



2012). In such situations, the energy extracted by the turbines is entrained mostly from 



the flow above the farm by turbulence, as quantified by Reynolds stress-mediated fluxes 



(Calaf et al. 


2010 


Cal et al. 


2010 



three-dimensional structure of these fluxes, in this section we investigate transport tubes 
of mass, momentum, and energy for eight different wind-turbine-array boundary layers, 
with different turbine spacings and configurations. 

Table [T] provides an overview of the different cases considered. Four cases use an aligned 
arrangement pattern, while four other cases use a staggered pattern (see Figure |4] for a 
sketch). These cases comprise different stream- wise spacings SxD and span- wise spacings 
SyD between turbines (with D the rotor diameter), as further detailed in Table 1. All 
results are obtained using LES following the approach discussed by 
Meyers fc Meneveau (2010). The four aligned cases are taken from 



Calaf et al. 



Calaf et al. 



(2010); 



( |2010[ ); 

the staggered cases are added in the current work, and are geometrically constructed by 
shifting every second span-wise row of turbines of the respective aligned cases along the 
span- wise direction. As a result, the stream- wise spacing between turbines doubles, while 
the span- wise spacing between rows is divided by two (see Figure |4]). 

The effect of wind-turbines in the LES is represented using an Actuator Disk Model 
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Figure 4: Aligned (left) and staggered (right) turbine arrangement patterns, and def- 
initions of stream-wise turbine-spacing Sx and span-wise turbine spacing Sy (non- 
dimensionalized by rotor diameter D). The mean flow [U) is in the x-direction 

(ADM). We consider cases with and without wake rotation. Cases without wake rotation 
do not include applied tangential forces at the turbine (ADM, i.e Case 1-8 and IF in Ta- 
ble 1). In another case (ADMR - Case IR in Table 1) we add tangential forces following 
the formulation used in Meyers & Meneveau ( 2010[ ). In a recent detailed validation study 
by |Wu fc Porte- Agel| pOl]| 7 it was demonstrated that except for near-wake effects close 
to the turbines with x < 3D, the non-rotating model (ADM) allows an accurate rep- 
resentation of the overall wake structures behind turbines. The rotating case (ADMR), 
including tangential forces, further improves near- wake behavior ( Wu fc Porte- Agel|2011 ). 
Moreover, in the same study, the Reynolds stresses were found to be accurately predicted 
by both formulations, thus allowing an accurate representation of the interaction of the 
wind farms with the atmospheric boundary layer. A snapshot of a typical LES velocity 
field using the ADM method is provided in Figure [5] Fu rther details on the methodology 
and computational set-up may be found inlCalaf et al. ( 2010[ ) and [Mey ers fc Meneveau 
(2010), and are summarized in Appendix |Aj where the effects of LES resolution (Case 1 
versus Case IF), and wake rotation (Case 1 versus Case IR) are also documented and 
discussed. 

4.1. Results: mean-flow stream tubes 

As a first step, we visualize classical stream tubes for Case 1. To this end, the three- 
dimensional velocity field is averaged in time to obtain a spatially periodic mean-flow 
velocity field with period corresponding to the turbine spacing, i.e. 7.85Z) in stream- wise, 
and 5.23D in span- wise directions for Case 1. Stream-tubes are obtained by constructing 
streamlines through 60 equally spaced seed points along a circle that coincides with the 
target turbine disk. The streamlines are tangent to the mean- velocity vector field obtained 
from LES (or later to the vector fields given by 2.2 or 2.9). During the procedure. 



we regularly add seed points whenever the curvature of downstream or upstream cross 
sections becomes too large, or stream lines are too widely spaced along the tube mantle. 

In Figure [6]ja) we show the mean stream-wise velocity field in a 8 by 1 turbine row, 
together with the turbine-rotor stream tube through the downstream turbine, i.e., defined 
by the streamlines through the rotor disk of that turbine. Note that the velocity field 
is periodic, but the stream tube is not. In Figure [6j^b) the energy tube is also shown 
for Case 1, illustrating large differences between both types of tubes. The differences 
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Table 1: Turbine arrangement parameters, and spacing of the cases considered for evalu- 
ating transport tubes in the current work (see Figure |4] for definition of stream- wise and 
span- wise turbine spacings). Number of turbines Nt in the simulation domain, domain 
size Lx ^ ^ 

H, and resolution N^x NyX of the LES are also shown (further details 
on the set-up, etc., are found in Calaf et al||2010 and in Appendix B). 











Sy 


Nt 


Lx ^ Ly X H 




NxXNyX 


turbine model 


Case 1 


alig 


ned 


7.85 


5.24 


8x6 


27r X TT X 1 




128x192x61 


ADM 


Case 2 


alig 


ned 


6.41 


6.41 


10x5 


2.047rxl.027r 


xl 


128x192x60 


ADM 


Case 3 


alig 


ned 


9.07 


4.54 


7x7 


2.027rxl.0l7r 


xl 


128x192x61 


ADM 


Case 4 


alig 


ned 


15.7 


10.5 


4x3 


27r X TT X 1 




128x192x61 


ADM 


Case 5 


sta 


ggered 


15.7 


2.62 


8x6 


27r X TT X 1 




128x192x61 


ADM 


Case 6 


sta 


ggered 


12.8 


3.21 


10x5 


2.047rxl.027r 


xl 


128x192x60 


ADM 


Case 7 


sta 


ggered 


18.1 


2.27 


8x7 


2.3l7rxl.0l7r 


xl 


128x192x61 


ADM 


Case 8 


sta 


ggered 


31.4 


5.24 


4x3 


27r X TT X 1 




128x192x61 


ADM 


Case IF 


alig 


ned 


7.85 


5.24 


8x6 


27r X TT X 1 




192x320x102 


ADM 


Case IR 


alig 


ned 


7.85 


5.24 


8x6 


27r X TT X 1 




128x192x61 


ADMR 




10 20 

Streamwise velocit}' 



Figure 5: Snapshot of a stream- wise velocity field in LES of a wind-turbine-array bound- 
ary layer. The color scale is stream-wise velocity in units of = [— {H / p)dpao / dx]^ ^"^ . 



are due to the considerable transport across the stream tube associated with turbulence 
(Reynolds stresses). Lebron et al. (20121 measured such fluxes through the mantle of 
a stream-tube using wind-tunnel data from a model wind farm, and found that the 
turbulent fluxes were dominant. Further discussion of this case is continued below, but 
first, we present stream tubes for the other wind- farm cases introduced in Table [T] To this 
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Figure 6: Mean stream-wise velocity field, (left) stream tube, and (right) total mechanical 
energy tube in a turbine row of a fully developed wind-turbine-array boundary layer flow 
(corresponding to Case 1 in Table [T]). The color scale is stream- wise velocity in units of 

= [-{H/p)dp^/dxY/^. 
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Figure 7: Upstream and downstream sections of stream tubes for wind farms with dif- 
ferent turbine spacings. (a-b) aligned Cases 1, and 4; and (c-d) staggered Cases 5 and 

8 - see Table [l] for details. ( ): turbine rotor; ( — ): sections at different upstream and 

downstream locations, with distances corresponding to a; = ±nSxD, and n = 2, 4, • • • ,10 
(farthest sections, aX n = ±10, are labeled). (•••): in (c)-(d) corresponds to the location 
of the staggered row of turbines. 



end, we display sections of rotor-disk stream tubes at upstream and downstream rotor 
planes in Figure [7]ja)-(d) for a selected number of wind-farm cases. We observe, again, 
that the stream tube continuously deforms further and further away from the rotor disk, 
and that the tube center does not remain at hub height. It is appreciated from these 
periodic cuts ( "Poincare sections" ) that the average mass flux through the turbine rotor 
plane is originating upstream from below the turbine level, while downstream it is ejected 
above turbine level. The main difference between the aligned cases (Figure[7]Ja),(b)) and 
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S 0.9 




Figure 8: (a,b) Evaluation of stream tubes in a wind- farm boundary layer with turbine 

spacing Sx — 7.85, and Sy = 5.23 with ( — ): Case 1, and ( ): Case IR (see Table 1). 

(a) Surface area along a turbine-rotor stream tube, (b) Flux of mean axial momentum 
through the stream tube, (c) Geometrical center of gravity of the stream tubes (all cases, 
see Table 1), with closed symbols: aligned case; and open symbols: staggered case. ■,□: 
Case 1, 5; o,»: Case 2, 6; ►,>: Case 3, 7; T, V: Case 4, 8; and x: Case IR. 



the staggered cases ([7](c),(d)), is that the stream tubes in the latter cases extend much less 
to the sides, as the sideways development of these tubes in constrained by the neighboring 
out-of-plane turbine rows. For sake of brevity, stream tubes of other cases are not shown 
here, as they have features which are very similar to the cases shown in Figure [7] 

Before turning to the main topic of momentum and energy transport tubes, first we 
report further properties of the classical stream tubes by evaluating the area {S, inversely 
proportional to the section-averaged mean velocity) of vertical cross sections of the tubes, 
as well as the axial fluxes crossing these sections as function of downstream distance. 
The evolution of S and the axial momentum flux are plotted for Case 1 and Case IR 
in Figure |8ja) and (b). It is observed that the evolution of cross-sectional surface and 
momentum flux along the tubes is quite similar for both cases, with and without rotation. 
Overall, we find that effects of wake rotation do not dominate transport of momentum or 
energy, and further discussion, comparing Case 1 and Case IR, is provided in AppendixjA] 

In Figure |8]J a) the surface 5' is displayed as function of the upstream and downstream 
distance from the tube's originating turbine disk. At a; = 0, it is observed that the slope 
of S is positive, associated with a reduction of the flow velocity by the turbine disk 
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thrust forces. Further downstream < x/{sxD) < 1, the surface decreases again, related 
to a speed-up of the flow (wake recovery). Also in upstream and downstream turbine 
planes (—6 < x/{sxD) < 6), similar trends are observed: at the turbine planes the area S 
increases (slow down of the flow), in between turbines it decreases. It is further observed 
that the difference between maximum and minimum S decreases farther from a; = 0, as 
the intersection of the stream tube with the turbine rotors at upstream and downstream 
planes decreases. Sufficiently far upstream (or downstream), e.g., x/{sxD) ^ —6, the 
trends change. Now the tube area S shrinks at turbine planes (with speed up of the 
flow), and grows in between (slow down). Here, the tube is not intersecting anymore with 
the turbine rotor. Hence, the average flow speeds up at the rotor, i.e. flow is partially 
driven around the rotor, and the flow slows down in between rotor planes, i.e. part of 
its momentum is transferred to the wake regions behind the turbine rotors by Reynolds- 
stress interactions. 

In Figure |8j(b) the flux of axial momentum through the stream tube is shown. Trends 
observed can be explained using the same rationale as above, and are largely related to 
the effective intersection between the stream tube and the turbine rotor disk regions. We 
further observe that the maxima of fluxes of momentum through the tube increase for 
—6 < x/{sxD) < 10. This is explained by the ascending trajectory of the stream tube, 
and the increase of mean-flow background momentum which is available to replenish 
momentum in the turbine wakes. Fluxes of energy through the stream tube, look very 
similar to the evolution of momentum fluxes, and are not further shown here. 

From the analysis above and Figure [8]ja,b), it is appreciated that the flux of axial mo- 
mentum through conventional stream tubes is highly non-trivial, affected by the upward 
motion of the tubes through the farm, together with Reynolds-stress exchanges over the 
tube mantle. Trends for the other cases (not shown) are the same. For all cases we flnd 
that through the mean velocity field, fluid volume (or mass) comes from below the tur- 
bines and downstream is ejected above the turbines. This is illustrated in Figure |8](c) , 
where the geometric center of the stream tube vertical cross sections is presented for all 
cases. We also observe that the geometric center for staggered cases remains closer to 
the turbine center. Note that the stream-wise turbine spacing Sx in the staggered cases 
is twice that of the aligned cases, so that this difference between staggered and aligned 
effect would be even more pronounced when the stream-wise distance is not normalized 
by SxD. 



4.2. Results: Momentum and energy tubes 

We now turn to the determination of momentum and energy transport tubes as defined 
in Sj2] Since the molecular viscosity in the LES is set to zero, and the contribution of the 
sub-grid eddy-viscosity compared to the resolved stresses is negligible (except very close 
to the ground), we consider only the Reynolds stresses based on the resolved velocity 
fleld in deflning the transport velocities. 

In Figures [9] and 10 sections of transport tubes of axial momentum and mean-flow 
mechanical energy, respectively, are shown for the various cases of Table I. These trans- 
port tubes are considerably different from the conventional stream tubes shown before. 
Neither momentum nor energy are conserved in the tubes, i.e. large sinks exist when 
the tubes (partially) pass an upstream turbine-rotor disk, such that the total tube cross- 
sectional area gradually shrinks until it reaches its originating turbine-rotor plane. At the 
originating rotor a large part of the remaining momentum/energy is removed, and due to 
effects of dissipation and further momentum/energy extraction at downstream turbines, 
the tubes rapidly shrink to zero. The additional one to two downstream tube sections 
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Figure 9: Upstream sections of axial-momentum transport tubes for wind farms with 
different turbine spacings. (a-h) corresponds respectively with Case 1-8, where (a-d) 

are all aligned cases, and (e-h) are staggered cases - see Table [T] for details. ( ): 

turbine rotor; ( — ): sections at different upstream locations, with upstream distances 
corresponding to x = —ns^D, and n = 2, 4, • • • ,20. (•••): in (e)-(h) corresponds to the 
location of the staggered row of turbines. 



that may be typically drawn before the tube disappears are not very enlightening; and 
therefore, only upstream sections are displayed for clarity. 

Furthermore, the energy tubes in Figure [TO] illustrate that the flux of total mechanical 
energy to the turbines in large wind-turbine-array boundary layers strongly depends 
on the stream- wise and span-wise spacing of turbines. In particular, when the span- 
wise spacing is large, energy is entrained from the sides, and only significantly further 
upstream is it entrained from above the turbines (see Fig. 10 i). In this case, the large 
span-wise spacing allows for high-speed flow to enter in between the turbine rows, where 
it further interacts sideways with the wake regions behind the turbines. For turbine arrays 
with narrow span- wise spacing (e.g. Fig. 10;, and all staggered cases. Fig. 10 3~h), it is 
observed that the energy is entrained directly from the flow above, and less from the 
sides. 
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Figure 10: Upstream sections of mean-flow mechanical energy transport tubes for wind 
farms with different turbine spacings. (a-h) corresponds respectively with Case 1-8, 
where (a-d) are all aligned cases, and (e-h) are staggered cases - see Table [T] for details. 
( ): turbine rotor; ( — ): sections at different upstream locations, with upstream dis- 
tances corresponding to a; = —nSxD, and n = 2, 4, • • • , 20. (• • • ): in (e)~(h) corresponds 
to the location of the staggered row of turbines. 



By investigating the energy flux in the tube in Figure 11 (here for Case 1), it is seen 
that the energy level remains largely constant between turbine planes (apart from effects 
of production of turbulence, which appears to be significant only in regions immediately 
downstream of the turbines, and power inserted by the driving pressure gradient), but 
drops significantly at the rotor disk locations as a result of the energy extraction by the 
turbines. This is quite different from the behavior of the stream tubes shown in Figure [8] 
We further find that results for Case 1 and Case IR (not shown) are essentially the same 
(cf. also Appendix [A|) . 

Finally, in Figure [l2[ a), the evolution of the energy-tube area is evaluated as function 
of upstream distance for the different wind-turbine arrays, while in Figure 12 b) the tube 
geometric center is displayed. Some differences between aligned and staggered cases are 
observed. In particular, the geometric center of the tubes moves upward more appreciably 
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Figure 11: Analysis of mean- flow mechanical energy tubes for Case 1: ( — ) flux of total 

energy (normalized by target turbine value); ( ) cumulative mean-flow dissipation 

by production of turbulence Pd{x); and (— •) cumulative power Poc{x) inserted by the 
driving force Vpoo- 
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Figure 12: (a) Evolution of cross-section area of energy tubes, and (b) geometric center 
Zcg of the energy tubes with closed symbols: aligned case; and open symbols: staggered 
case. ■,□: Case 1, 5; o,»: Case 2, 6; ►,>: Case 3, 7; and T, V: Case 4, 8. 



for the staggered cases - due to their narrower span-wise spacing, there is less room for 
sideways expansion. 



5. Discussion 

The visualization of mechanical-energy tubes, as shown in Figure |10[ provides an in- 
tuitive understanding about the region of the flow supplying the power that ultimately 
is available at individual turbines. In this section, further features of the energy tubes 
are discussed. 

First, it is useful to recall that the tube section associated with upstream plane 



x/{sxD) = —n (e.g., in Figure 10) also corresponds to the section in the plane a; = of 
a tube associated with a turbine at downstream distance x/{sxD) = n. This is a direct 
consequence of the periodic nature of the flux vector field. Now, recall that changes in 
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energy flux along the tubes are only due to internal sources and sinks, i.e. the power 
inserted by the driving force Vpoo (^oo)j the power extracted by turbines (Pt), a-nd the 
production of turbulence (-Pd). Hence, except for the relatively small difference between 
Pd and Poo compared to Pt, the sections in Figure 10 (omitting the minus signs in their 
numbering) , may be roughly interpreted as containing the total flux of mechanical energy 
in the current cross section of the wind farm that will be extracted away in the next n 
turbine rows. 

Furthermore, considering the subsequent sections of energy tubes in Figure |10[ we 
observe for all cases that the side and bottom boundaries of the sections asymptotically 
converge for increasing n. This suggests that the total boundary layer region is divided 
into a region covered by turbine tubes, and a region which is not. We investigate this in 
detail by seeding the cross section of the boundary layer with a large number of points, 
and constructing the downstream flux lines over a large number of periodic cycles. We first 
focus on Case 1 in Figure [TS] 100 by 100 points are seeded on a Cartesian grid covering 
a cross section of the boundary layer. Subsequently flux lines are tracked over a large 
number of cycles to determine their attractors. Making connections to dynamical systems, 
we remark that the flow corresponding to the energy flux lines is not conservative, and so 
unlike the velocity field which does not possess attractors due to volume conservation, the 



energy flux lines can have attractors and repellers. In Figure 13 [a,), the white (uncolored) 



area, represents all points from which the flux lines are attracted to the ground surface. 
The gray shaded area shows the points whose energy flux lines will eventually pass 
through a wind turbine disk. These flux lines are all attracted to what is seen as a point 
in the current cross (Poincare) section. To illustrate representative trajectories of flux 
lines from different areas, we selected four points of which the flux lines (projected onto 
the turbine plane) are drawn. From this it is seen that points from white areas are indeed 
attracted by the ground surface, while points from the gray area end up in a flux line 
close to the turbine center. 

In Figure [Isj^b) , we show the energy flux lines from the 2-D vector held obtained 
from steam- wise averaging the energy-flux vector fleld according to {Fe)x- This flgure 
illustrates, in an integrated sense, how the energy flux lines are attracted to different 
regions in the boundary layer. Note that the comparison with Fig. [iSj^a) is qualitative 
only, and not exact: the lines in Figure 13 ^a) are obtained through Lagrangian tracking. 
The effective flux field is changing in x over one period due to the expansion of the flow 
close to the turbines - these effects are however relatively small, such that a comparison 
remains insightful. 

Basins of attraction for most other cases are similar to Case 1. For instance, in Fig- 
'&) the basins of attraction for Case 3 are shown. As for Case 1, a relatively large 
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part of the flow total energy is "attracted" towards the turbine, while in-between tur- 
bine rows, a smaller area is attracted to the ground (but dissipated before reaching it 
since no work is being performed at the bottom boundary). We find one case where this 
picture differs significantly, i.e. Case 4 with Sx — 15.7 and Sy = 10.5, which is shown in 
Figure [Tgb). This is the case with the widest span-wise turbine spacing (and together 
with Case 8 has the highest average turbine spacing considered) . Three distinct domains 
are now observed: a domain where flux lines are attracted to the ground (in white), and 
two gray domains where the flux lines are attracted along spirals to the top left and top 
right of the turbine row. Finally, flux lines in the lighter gray domain flrst pass through 
a turbine disk, where energy is extracted. Flux lines in the darker gray domain do not 
pass through a turbine. The result in Figure [l4{^b) is interesting, as it suggests that the 
available driving power for this case is not maximally used for wind-energy conversion. 
In large domains of the boundary layer (i.e. the white and dark-gray domains) , the driv- 
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Figure 13: (a) Domains of attraction of energy-flux lines through a cross section of the 
boundary layer for Case 1. White domain: flux lines are attracted to the ground surface; 
gray domain: flux lines pass through a turbine, and are attracted to a line close to the 
turbine center. Full line: four different flux lines and their projected trajectory (starting 
from points marked with DjOjAjO respectively), (b) Projected flux lines in a turbine 
plane, obtained from the stream-wise averaged flux vector field {Fe)x (for Case 1). 

ing power is balanced by dissipation through production of turbulence, and only in the 
light-gray area is driving power converted to useful energy. 



6. Conclusions 

In the current work we explore the concept of momentum and energy transport tubes 
to study the three-dimensional mean fluxes of these properties in turbulent flows. These 
tubes are constructed based on transport vector fields, which include effects of Reynolds 
stresses and mean-flow viscous stresses. In particular, when transport processes are dom- 
inated by Reynolds stresses and turbulence instead of mean-flow convection, such trans- 
port tubes are an interesting means of visualizing where the momentum or the energy in 
the flow originate and/or are transported to. 

As illustrative application, we study stream-wise axial momentum and energy tubes 
in fully developed wind-turbine array boundary layers. Analyzing conventional stream 
tubes, we find that on average, the mean flow volume (mass) passing through a wind 
turbine disk comes from below the turbine, and is downstream ejected into layers above 
the turbines. Based on the energy tubes, we find that the energy takes a different path 
to reach the wind turbine locations. Depending on turbine arrangement, there are two 
distinct paths and mechanisms taken by the energy as it reaches the turbines: a sideways 
flux and a top-down flux. Sideways fluxes themselves are fed by a top-down flux in 
regions outside the turbine wake area. For large span-wise turbine spacings, sideways 
fluxes of energy dominate; for small span-wise spacings, only the top-down mechanism 
is dominant. 

Further investigating total-mechanical energy tubes and flux paths, we find that they 
define different basins of attraction in the boundary layer. In one part of the boundary 
layer, energy flux lines are attracted to the ground surface, while around and above 
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Figure 14: Domains (basins) of attraction of energy-flux lines through a cross section of 
the boundary layer, (a) Case 3 - White domain: flux lines are attracted to the ground 
surface; gray domain: flux lines pass through a turbine, and are attracted to a line close 
to the turbine center, (b) Case 4 - White domain: flux lines are attracted to the ground 
surface; gray domains: flux lines are attracted into a spiral, either on top left or top right 
of the turbine; light gray domain: flux lines pass though turbine disk; dark gray: flux lines 
do not pass through turbine disk. Full line: four different flux lines and their projected 
trajectory (starting from points marked with 0,0, A, respectively). 



turbines, flux lines are typically attracted to the turbine disk region. In some cases, 
attracting points above the wind turbine region were also observed. The relative size of 
these domains of attraction depend on turbine spacing and arrangement pattern. 

In closing, we point out that it also may be interesting to consider momentum and 
energy flux lines and tubes in the case of laminar flows, where only the viscous fluxes 
provide differences to regular velocity and stream-lines (some simple examples were pro- 
vided in Section |3]). Moreover, there exists the possibility of fruitful analogies between 
physical-space trajectories of the generalized transport vector field and trajectories in 
phase-space of dissipative dynamical systems, similarly to how Hamiltonian dynamical 



systems provide useful analogies for laminar-flow chaotic mixing ( Aref 1984 Ottino 1989). 
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Appendix A. Wind-farm simulations: governing equations and 
computational setup 

This appendix provides additional details about the LES. The methodology has already 
been presented in Calaf et al. (2010). Moreover Case 1-4 (see Table [I]) in the current 
work correspond with Cases A3, K, J, and G of Calaf et al. (2010), respectively. 

The current work visualizes the interactions between an infinite wind farm and a 
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neutral boundary layer. We simulate rough-wall fully-developed boundary-layers in a 
periodic domain, driven by a pressure gradient. We do not include Coriolis forces. While 
these are present in the atmosphere at the larger scales of the flow, the main rationale 
underlying this approach (see discussion in Calaf et al. ( 2010[ )) is based on the classical 
hypothesis that inner-layer dynamics of a boundary layer [y < 0.1577) are approximately 
independent of outer layer effects. We presume that turbines (with height sa 100m) are 
situated in the inner layer of the boundary layer, which is relevant in practice for many 
atmospheric cases with boundary layer depths above 1 km. This approach was used before 



by Calaf et al. (2010), where it allowed for the characterization of increased surface 



roughness induced by a wind farms, and the derivation of algebraic surface-roughness 
models. Later, wind-farm performance obtained based on these surface-roughness models 
were shown to compare well with observations in the Horns Rev, and Nysted wind farms 
( [Meyers fc Meneveau||2012[ ). 



A.l. Large-eddy simulations 

We consider thermally neutral flow that is driven by an imposed pressure gradient. LES 
solves the the filtered incompressible Navier-Stokes equations for neutral flows and the 
continuity equation, i.e.. 



duj 

Ik 



dui 
dxi 
dujUj 
dx-j 







1 dp ^ dr., 
P dx^ 



dx.. 



(Al) 
(A 2) 



where Ui is the resolved velocity field, p the pressure, Tij are the subgrid-scale stresses, 
and where the density p is assumed to remain constant. Furthermore, fi represents forces 
introduced by the turbines on the flow (see discussion below). Since the Reynolds number 
in atmospheric boundary layers away from the bottom boundary is very high, we neglect 
the resolved effects of viscous stresses in the LES. The deviatoric part of the subgrid-scale 
stresses is modeled here with the conventional Smagorinsky ( 1963| ) model, with a constant 
coefficient Cs=0.14 (the trace of the subgrid-scale stresses Tfcfe/S is not modeled, but 
instead absorbed into the pressure term, as is common practice in LES of incompressible 
fiow). Near the bottom surface, the Smagorinsky length scale A (= C^A far from the 



surface) is damped using the classic wall damping function of Mason & Thomson ( 1992 ), 
i.e. A~" = [CsA]~" -I- [k{z + ^;o,;o)]~", where we take n = 3. Other works have also 
used more advanced subgrid-scale models in LES of wind farms (e.g. the scale-dependent 
Lagrangian model of |Bou-Zeid, Meneveau fc Parlange||2005| was used for several of the 
simulations presented in Calaf et al. (20101), but the differences in mean velocity and 
Reynolds stress distributions were found small, especially in regions where the transport 
tubes of interest in this study are mostly located. 

In the stream-wise and span-wise directions, we use periodic boundary conditions. 
The top boundary uses zero vertical velocity and zero shear stress condition. At the 
bottom surface, we impose zero normal velocity and use a classic, imposed wall-stress 
boundary condition. It relates the wall stress to the velocity at the first grid-point using 
the standard log (Monin-Obukhov) similarity law ( Moeng|[l984 ): 



\n{z/zQ^io) 



0.5 , 



Tw2 



(A3) 



(A 4) 
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where the hat on u and v represents a local average obtained by filtering the LES velocity 
field with filter width 4A (see |Bou-Zeid et a^.| ( ,2005J for more details about such filtering). 

The simulation code uses a pseudo-spectral discretization in the horizontal directions. 
The nonlinear convective terms and the SGS stress are de-aliased using the 3/2 rule 
( Canuto et a/.||198"8 |. Message Passing Interface (MPI) is used to run the simulations in 



parallel mode, and the FFTW library is employed for Fourier transforms ( Frigo & John- 



son 



2005). In the vertical direction, a fourth-order energy-conservative finite-difference 
discretization is used ( Verstappen fc Veldman|2003 ) . Time-integration is performed using 
a classical four-stage fourth-order Runge-Kutta scheme. 



A. 2. Turbine forces 

For the wind-turbine forces, we use an actuator disk model. These type of turbine rep- 



resentations were adopted in LES by Jimenez et al. (2007 2008); Ivanell et al. (2009); 



Calaf et al. (2010) amongst others. Recently, they were thoroughly vahdated against 



wind-tunnel data by Wu & Porte- Agel (2011 1. Details of the current implementation are 
given in Meyers & Meneveau (2010). In this model, the thrust and tangential forces in 



the turbine rotor disk per unit actuator-disk area are given by 



Ft = 



Crp 



Fe(r) 



1 w ,-T\2 {u)d 



(A 5) 
(A 6) 



with the subscript 'd' denoting an averaging over the turbine disk region, and the super- 
script 'T' denoting time filtering or averaging over a time-scale of order T. Thus {u^)d 
is the disk averaged and time-filtered velocity (further discussed below). The parameters 
C'rp and C'p are modified thrust and power coefficients, defined based on the turbine disk 
velocity instead of the undisturbed upstream velocity as conventionally used. Their val- 
ues are directly related to the aerodynamic lift and drag coefficients of the turbine blades, 
the blade geometry, etc. ( Meyers fc Meneveau|2010 1. We use values of = 4/3, and for 
the rotating (ADMR) case C'p = 1 (for the non-rotating case, the tangential forces are 
zero). For a lone-standing turbine, these values would correspond to conventional thrust 
and power coefficients of Ct = 0.75 and Cp w 0.42 ( [Meyers fc Meneveau|[20Tol ) . Finally, 
n is the turbine angular velocity, and < r < D /2 the radial location on the turbine 
disk. For the ADMR case in the current work (CaselR), we use flR/u^, = 60, which 
roughly corresponds to a tip-speed ratio of A w 6.7 (based on the average boundary-layer 
velocity at turbine hub height). 



To implement the forces fi in (A2), the turbine forces Fi (with axial, and tangential 
components Ft, and Fg) are first described in the turbine- rotor plane. In a second step, 
these forces are filtered using a Gaussian convolution filter on locations which corre- 
spond with the coordinates of the LES grid. We use a Gaussian filter, with filter width 
A = 1.5A (and A is the grid spacing) to avoid Gibbs oscillations on the LES grid. A sim- 
ilar smoothing approach was used for actuator-line representations in |Sorensen fc Shen] 
(2002). To evaluate the disk averaged local velocity {u)d needed for the determination 
of the force Ft and Fg, we employ the geometrical rotor footprint on the LES grid as a 
weighting function for the averaging. Moreover, is obtained from ({t)^ by using a 

one-sided exponential time-filter, using a time window of Tu^,/H = 0.6. Further details 



are found in Calaf et al. (2010) and Meyers fc Meneveau (2010). 
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Figure 15: (a) Upstream and downstream sections of stream tubes (at x = zL^s^D and 
±10sxD) and (b) upstream sections of total mechanical energy tubes (at x = —nSxD 
with 71 = 5, 10, 15, 20) for (— ): Case 1; and (— ): Case IF — see Table[l|for details. 
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Figure 16: Upstream and downstream sections of stream tubes for (a) non-rotating case 

(Case 1); and (b) rotating case (Case IR) — see Tablejljfor details. ( ): turbine rotor; 

( — ): sections at a; = iSsj^D; and (— •): sections ai x — ±.\QsxD. 



A. 3. Effects of numerical discretization 

All cases except Case IF are discretized using similar domain and grid sizes. For the 
steam-wise and span-wise domain size we use — 2ttH and Ly ~ ttH (with H the 
boundary-layer depth and domain height) to allow for the large scale turbulent structures 
that typically emerge in the boundary layers to be represented properly. The cell size in 
the simulations correspond to Ay « Az « 0.016i/, and Ax w 0.05H, except for Case IF, 
where a finer mesh is used, with Ay w Az w O.OIH, and Ax « 0.033-ff. The effect of 
refining the mesh on stream tubes and mechanical energy tubes is shown in Figure [TSl We 
observe that the shape of the tubes remains relatively unaffected, especially considering 
that the tube cross-sections shown correspond to very long integration distances of the 
stream and energy lines and small differences in mean velocity arising from the LES can 
accumulate while integrating the trajectories. 

A. 4. Effects of wake rotation 



In Figure 16 we examine effects of wake angular momentum introduced via tangential 
forces on the evolution of the stream tubes, by comparing Case 1 and Case IR. As 
expected, the induced swirl is associated with a twist of the stream tubes in upstream and 
downstream directions when compared to the symmetric non-rotating case. Nevertheless, 
also in the rotating case, mass is ejected upwards downstream from the turbines, while 
upstream it is entrained from below the turbines. Moreover, as already briefly discussed 
in Section |4.H momentum fluxes and cross sectional area along the tube are not much 
influenced by wake rotation. 

In Figure [17] the momentum and total mechanical energy tubes are displayed for Case 
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Figure 17: Upstream sections of (a) axial momentum, and (b) mean-flow mechanical en- 
ergy transport tubes for Case IR. ( ): turbine rotor; ( — ): sections at diiTerent upstream 

locations, with upstream distances corresponding io x — —nSxD, and n = 2, 4, • • • ,20. 



IR. When comparing with the non-rotating cases, i.e. Figure |9]ja) momentum, and Fig- 
ure [Tola) for mean-flow mechanical energy, it is appreciated that the effect of rotation 



appears less visible than for the conventional stream tubes (cf. Figure 16). In particular, 
the effect of rotation is quite small for energy tubes. This is not unexpected, as modern 
wind farms operate at relatively high tip-speed ratios, leading to low torque (for a given 
amount of power) and low associated tangential forces, such that effects on the Reynolds 
stresses, which are largely responsible for transport of energy, remain small. 
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